The perturbative scalar massless propagator in Schwarzschild spacetime 
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A short derivation is given of the weak gravitational field approximation to the scalar massless 
propagator in Schwarzschild spacetime obtained by Paszko using the path-integral approach. The 
contribution from the direct coupling of the quantum field to the scalar curvature is explictly in- 
cluded. The propagator complies with Hadamard's pattern, and the vacuum state is consistent with 
the perturbative version of the Boulware vacuum. The momentum space propagator is computed 
for massless or massive particles to the same perturbative order. The renormalized value of {(jP'ix)) 
for the massless case is reproduced. 

PACS numbers: 04.62. +v 

In [l[ , Paszko carried out a perturbative calculation of the Feynman propagator of a scalar particle in presence of 
a Schwarzschild metric, to first order in the black hole mass, obtaining an explicit analytical expression for massless 
particles. A path- integral approach was adopted for that calculation. In this note we show that this result can be 
obtained with considerably less effort by means of the operator formalism. 

For simplicity we will work out the Euclidean version and perform the Wick rotation at the end to obtain the 
Feynman propagator. We use units h = c = I and Af , with dimension of length, denotes the black hole mass times 
Newton's constant. 

Our starting equation is 

{-a + m^ +^R)^G{x,x') }. 8{x-x'). (1) 

V5(2;) 

Here □ = g~^l'^d^g^l^g^^ is the Laplacian operator, R the scalar curvature, and G{x, x') is the (Euclidean) 
propagator from x' to x of the scalar particle with mass m. 

Following [l|, for the (Euclidean) Schwarzschild metric with mass M, we use isotropic coordinates [1], 
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Here we use r to denote the Euclidean time (not the proper time) , in order to distinguish it from the Lorentzian time 
t, below. 

To first order in an expansion in powers of M the propagator equation reduces to 
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G{x, a-') = ( 1 - ^ j S{x - x') + 0{Ar), 



(3) 



Let us expand also the propagator and the scalar curvature, 

G = G("' + G(i' + 0(m2), R = + 0{M^), (4) 

where G^") and contain precisely n powers of M . Substituting in Eq. ([3]) and equating terms with equal powers 
of M produces the following relations 

(-a^- V2 + m2)G(°)(x,x') = S{x~x'), (5) 



i^d^-V^ + m^)G^^\x,x') 
From these relations it immediately follows that^ 

G^^\x,x') = I d^x"G^''\x,x") 



AM 



4Af 



{dl„^\m^)-iR^^\x") 



G(")(a;,a:'). 
G(")(^",a;')- 



(6) 



(7) 



Once the Minkowski vacuum has been adopted at zeroth order, the perturbative expansion automatically selects a privileged choice of 
boundary conditions at any other perturbative order. We do not override that choice here. 
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FIG. 1: Feynman graph for Eq. ((7|. The vertex can be read ofT from Eq. ([8]). 

Of course, the same result can be obtained by applying standard perturbation theory to first order. See Fig[TJ The 
Feynman rule for the interaction vertex of the particle with the static gravitational field can be read off from the 
Lagrangian: 

= \{d,c^f + \nr'<t>' + \^ {{dr^f + ^-^0^) + l^R^'U' + O(Af^). (8) 

Note that in the Euclidean perturbative computation the metric is just 6^^,. 

In order to obtain a closed expression in x-space (as opposed to momentum space), from now on we consider only 
the massless case, m — Q. In this case 

As follows from Eq. ([7]), the effect of ^i? on the propagator is additive to first order, so we define 

G^^\x,x')^Gi^\x,x')+Gf\x,x'). (10) 

Let us consider first the effect of the scalar curvature term. From Einstein's equations, a Dirac delta mass distribu- 
tion as source of the Schwarzschild metric corresponds to a scalar curvature R{x) = 8ttMS{x) + 0{AP). When this 
is inserted in Eq. (O, one readily obtains 



Gt (a;, X ) — C \9 / A 9 I / I /\9 ' (-'--'-) 
^ (2n)'' rr i\T'^ + (r + r'Y 

where Ar — t — t' . For the ^ = part, Eq. ([7]) reduces to 

g[,^' {x, x') = AMdlF{x, x'), (12) 



with 



{2^fF{x,x') = / d^x"\ ^ „ ^ (13) 

r" (x — X Y (x — x Y 



Here it is appropriate to use a standard Feynman parameterization 
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ab 

which yields 



— = / dC 2' (14) 
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where y = (x + (1 — C)x' and P = (a; — x')"^. Now it is straightforward to carry out the x" integration, starting with 
the angular part, then r" and lastly r" . This gives 



Jo \y\ \Vy' + Cii-Oi'-\y\) 

The remaining integral produces a rather complicated expression. However, the function F itself is not needed. As it 
turns out, the first time derivative of F takes a simple form: 

{4TT)^drF{x,x') = dC , ^ 

Jo vc^' + (1 - cy^ + c(i - c)At2 



4 ^ At 



arctan ( ) . (17) 



Taking another time derivative yields the final expression for G*-^' (adding the ^ dependent term): 
-,(1), ,x 1 4M /2At , f At \ r + r' f , „ l'^ 



G};'{x,x') = jT-^-pT -p^arctsin - -— 1 + ^— . (18) 

^ ^ ' ' (27r)2 P \ P \r + r'J (At)^ + (r + r')^ \ 2rr' J J ^ 

To obtain the propagator with Lorcntzian signature, it only remains to perform the Wick rotation. This entails 
taking an analytical continuation 

iApiAt- X, x') = G(Ar = iAt; x, x'). (19) 

Specifically, for the Feynman propagator, for positive (negative) At = t — t' , the continuation is to be taken from the 
positive (negative) Re(Ar) half-planes, respectively. Hence, At = iAt + sgn{At)ri, rj ~ 0+. For the zeroth order, this 
yields the free Feynman propagator for a massless scalar particle 

aP{x,x')= ' / . , (20) 

where = At^ — Ax^ with Ax = x — x' . For the first order term, the Wick rotation gives (leaving 77 implicit) 
(1)^ i AM / r + r' ( , ^ \ IAt , ( At 



A^j;'{x,x') = 77-^^ -r-5 7 7W 1 - Ctt-T + ^arctanh . (21) 

' ' (27r)2 s2 yA<2-(r + r')2 \ ^ 2rr' J \r + r' J J ^ ' 

This form applies to the case |At| < r + r' (taking for arctanh the real branch). When | At| > r + r', the expression 

is2 

„ i AM { r + r' ( \ 2At , fr + r'\ \At\\ , , 

A^p\x,x') = , .„ „ — 1 - i— + ^arctanh — — - ivr^ . 22 

^ (27r)2 s2 \At^-(r + r'Y \ 2rr' ) \ At J J 



When ^ = 0, this is the same result already obtained in [l|. 

i'-°\x,x') + A'-^ 



We have explictly verified that Ap\x,x') + Ap\x,x') solves the Klein-Gordon differential equation to 0{M). 



From the symmetries of the problem, the propagator will be a function of r and r', as well as and At. To 0{M), 
the ^ = part of the propagator depends only on the combination r -I- r' , but this property is not preserved by the 
contribution from ^i?. 

To the same order in M, the explicit result in momentum space can also be given, even for massive particles, as it 
follows immediately from the Feynman vertex in Eq. ([8]) (in its Lorentzian signature version). Letting 

D{p) = ^-T--^ (23) 



This is most easily obtained by using first the identity arctan (i) = — arctan(x) -|- ■^sgn(x) in Eq. JTs}, and noting that under Wick 
rotation sgn(AT) goes to sgn(At). 
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the propagator is given by 



{2t:)H{p-p')D{p), 
2t:S{p° ~p'°)8nMD{p) 



2p°p'° 



Next, we investigate whether the propagator conforms to Hadaniard's pattern 



D{p'). 



2i{2TTfl\F{x,x ) ^ - + yiog(cr) 

(7 



(24) 



(25) 



where J7, V and W are regular functions of x and a;' at a; = x' (at the points where the metric itseh is non singular), 
with U{x,x) = 1, and a{x,x') is Synge's function, half the square of the length along the geodesic joining x and x' 
(and cr —> cr — z?7 is implicit). To be a Hadamard state is a well-known requirement for the vacuum being sufficiently 
regular in the ultraviolet region to ensure a non singular stress-energy tensor [3, IH . 

The momenta q ~ p — p' and Q ~ {p + p')/2 control, respectively, (x + x')/2 and x — x' in x-spacc. The large Q 



behavior A^-* ~ Q'^^/Q^ competes with that of A^'' ~ l/Q^, and this could give rise to a singular coincidence limit 
in the propagator in x-space. The corresponding situation manifests in s'^A^p\x,x'), which behaves as At^ / and 
so it is not a continuous function in the coincidence limit. 

In order to clarify this issue, a has to be computed to order M . Being an extremal curve, the length of the geodesic 
can be obtained to first order by using the zeroth order path, i.e., the straight line. This gives 



,(0) 



a{x,x') = is2-A/(Ai2 



where h is the average value of 1 /r along the path; 



h{x^ x') ~ 

Jo 



dX 



Aa; + (1 - A)a;'| |Aa; 



■arctanh 



\Ax[ 

r + r' 



(26) 



(27) 



When combined with the propagator in Eq. (|2T|) . the following relation obtains 



i{2TTf2aAF = 1 + 8A/ 



Ai 



2 / arctanh ( | arctanh ( 



-h4M 



Ai2 - (t 



2rr' 



\Ax\ 
arctanh 




0(M2). 



(28) 



The term with 4M is manifestly regular as a function of Ax'' at x = x' . In the term with 8M, the structure 
(/(At^) — f{Ax^))/{At^ — Aa;^), being / a regular function, is itself regular. It follows that the propagator is of 
Hadamard type with V = 0{M^), and U{x,x) = I + 0{RP) is also checked. 

To further expose the vacuum structure, we adopt the point of view that the Minkowski vacuum, |0), and the 
corresponding Hamiltonian Hq, arc subject to an 0{M) time- independent perturbation. Specifically, H = Hq + Hi + 
0{A'P) with (the canonical momentum being 7r(a:) = (1 + 4,M / r)dQ(j){x) + 0{AP)) 



Hi = / d?x 



1 



4Af , 1 



m2(/)2(a;) - tt'^{x)) + £,R'^^\x)4)'^ {x) 



(29) 



Straightforward application of timc-indcpcndcnt perturbation theory to first order yields the new perturbed ground 
state: 



|0) = |0) 



SttM 



{27r)3 V2^ {27t)3 V2^ uj + uj' 



Ife-hfc'P 



a){k)a}{k')\Q)+0{AP), 



(30) 



where 



and I 



^, and at(fc) is the creation operator of the unperturbed Fock vacuum, with 



normalization [a{k) , {k')] = (27r)^(5(fc - k'). 

The perturbed vacuum is annihilated by the perturbed annihilation operators 



a{k) ~ a{k) 



d^k' 
(2^ 



1 



1 2 
2™ 



- OJUJ 



a{k') 



SttAI 

UJ + Uj' 



\k + k' 



a\k')] +0(A/2). (31) 
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These operators satisfy [H, d{k)] = —u>a{k) to 0{M). The single particle spectrum is unchanged since it is determined 
by the asymptotic scattering region, far from the black hole. The fact that a{k)\0) = and the commutator of the 
perturbed creation and annihilation operators are still c-numbers allows to use Wick's theorem for the field correlation 
functions, and this implies that the state |0) is Gaussian. 

A direct calculation shows that the Fourier transform of — z(O|T((/)(x)0(a;'))|O) just reproduces the propagator in 
Eq. ([24]) . (Note that H, and not Hq, has to be used here to evolve the fields in Heisenberg image.) The previous 
arguments strongly suggest that |0)is just the Boulware vacuum, the stationary non thermal and non radiating ground 
state of the Schwarzschild metric [a, 01: albeit to 0{M). Therefore, the propagator applies to a static spherical star, 
in the regime r ^ M, rather than to a proper black hole, which would be in a thermal state with a temperature 
T = l/(87rM) 0. 

The renormalized value of {(p{x)) has been computed in Q for a massless particle in the weak gravitational field 
of a static spherical star. Local approximations have been avoided in Q. 

From the definition of the propagator, the bare value of {(fP{x)) follows from taking the coincidence limit in G{x, x') 
or iApix^x'). Here we take the alternative approach of going back to Eq. ([7]) and setting x ^ x' there (for m = 0). 
Straightforward integration over r" gives 



G^^\x,x) 



irM 



d^x" 



1 



1 



\x — X 



//|5 



1 



\x — X 



lf\3 



(32) 



The term with ^ is ultraviolet finite and it just reproduces the result of taking the coincidence limit directly in Eq. (jlip . 
The term without ^ is divergent. However, as a distribution \x — x"\~^ — (l/6)V^|a; — x"\^^. Integration by parts 



and use of V la;'' 



'4:TtS{x") immediately gives 



{^\x)) 



1 M 



0{M^). 



(33) 



This result is in agreement with that obtained in [9| by using dimensional regularization plus standard subtraction of 
the singular part [2| . 

The Feynman propagator at x = x' but arbitrary At, for the same setting (weak static spherical gravitational field 
and massless quantum field) follows from the results in (lo| . Use of i?(i)(a;) = 8ttMS{x) in Eq. (6.9) of [H, yields 
a result that almost agrees with our value for A'-p\t,x]t' ,x) from Eq. The resuh in ^ is non singular in 

the limit At — >■ 0, and it is consistent with Eq. (j33p . Instead, we obtain the same non singular part plus a divergent 
contribution: 



iA'p'(t,x;t',x) 



div 



1 2M 1 



AP' 



(34) 



This expression is consistent with Eq. (Al) of [Tll |. 

As a final comment, we note that ultraviolet divergences are expected to arise at higher orders in a strict expansion of 
the propagator in powers of M . These would come from small radii in the intermediate point integrations. Inspection 
of the integrals involved suggests that such an expansion breaks down already at 0{M^) due to the presence of terms 
{M/r"f. 
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